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We have derived a variational principle that defines the nonequilibrium steady state transport 
across a correlated impurity (mimicking e.g. a quantum dot) coupled to biased leads. This varia- 
tional principle has been specialized to a Gutzwiller's variational space, and applied to the study of 
the simple single-orbital Anderson impurity model at half-filling, finding a good qualitative accord 
with the observed behavior in quantum dots for the expected regime of values of the bias. Beyond 
the purely theoretical interest in the formal definition of a variational principle in a nonequilibrium 
problem, the particular methods proposed have the important advantage to be simple and flexible 
enough to deal with more complicated systems and variational spaces. 
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I. INTRODUCTION 

Nanocontacts of quantum dots, single molecules or 
atoms, and nanowires are ideal candidates to realize elec- 
tronic devices where a source-drain current across the 
contact can be magnetically controlled. Indeed, because 
of the low dimensionality of the contact region, elec- 
tronic correlations grow in strength and may stabilize a 
local magnetism that influences electron tunneling. The 
Kondo-like zero-bias anomalies first observed in quantum 
dotsi are just the simplest manifestation of such a local 
magnetism, but one can foresee even more spectacular 
phenomena, like giant magnetoconductance^ 

From the theory side, this is a complicated problem 
first of all because electronic correlation is the main actor 
and is difficult to treat, and secondly because the inelastic 
tunneling spectrum requires full out-of-equilibrium calcu- 
lations. Many complementary techniques have been used 
to characterize the nanocontact at equilibrium. For in- 
stance ab initio Local Density Approximation (LDA) cal- 
culations can provide the electronic structure and predict 
whether magnetism could indeed be stabilized^— at least 
at the mean field level. Inclusion of quantum fluctuations 
requires many-body techniques, like numerical rcnormal- 
ization group^^ which are often applied to oversimplified 
models, like the single-orbital Anderson impurity model, 
although there are recent attempts to join together the 
two approaches Unfortunately, out of equilibrium 
properties are much more difficult to study. Apart from 
many-body Keldish perturbation theory,— many sophis- 
ticated numerical techniques have been developed in re- 
cent years to cope simultaneously with out-of-equilibrium 
and strong correlations*^— However, given the com- 
plexity of the electronic structure that may arise at a 
nanocontact e.g. of a molecule or a bridging transition 
metal atom, it would be desirable to have at disposal ap- 
proximate techniques simple and flexible enough to deal 
with realistic situations otherwise prohibitive with more 
accurate numerical approaches, as those previously men- 
tioned. 

In this paper we shall propose an out-of-equilibrium 



extension of the conventional variational approach, and, 
on such basis, an out-of-equilibrium extension of the 
Gutzwiller approximatio n 17 ' 18 for correlated electron sys- 
tems. 

The paper is organized as follows. In Sec. [II] we briefly 
introduce the Hershficld formulation of the nonequilib- 
rium steady state problem in quantum dots. In Sec. IIIII 
we introduce the concept of scattering operators and 
derive some results related with the resonant model. 
In Sec. IIVI and [V] we derive the variational principle 
that defines the nonequilibrium steady state transport 
across a finite junction (e.g. a quantum dot) coupled 
to biased infinite leads. In Sec. I VII we formulate a 
Fermi-liquid assumption for the system in the low en- 
crgy/temperature/bias regime. In Sec. IVIll we introduce 
very briefly the standard Gutzwiller variational method 
for the single band Anderson impurity model in equilib- 
rium. In Sec. IVIIII we propose a generalization of the 
Gutzwiller variational method to nonequilibrium. Fi- 
nally, Sec. IIXI is devoted to the conclusions. 



II. THE PROBLEM 

We consider two biased macroscopic leads described 
by non-interacting electrons coupled to a bridging region, 
the quantum dot, described by discrete electronic multi- 
plets 



n = T + v + n l 



(i) 



where V describes the tunnelling between the leads and 
the nanocontact and Hmt the local interaction in the 
nanocontact. 

One assumes that initially the leads are not coupled 
through the bridging region, each lead being subject to a 
different electrochemical potential. Such a situation can 
be described by a density matrix 



Po 



-|Stto(*) 



/Tr(e' 



where 



Ho($) = Ho + $Y , 



(2) 



(3) 



2 



with Ho the non-interacting Hamiltonian of the indepen- 
dent left (L) and right (R) leads plus the nanocontact 

Ho = f + n mt , (4) 

$ is the applied voltage between the two leads, and 



$Y = $(N L -N R )/2 



(5) 



that describes the electrostatic energy gain due to the 
presence of the bias voltage, where N is the number 
operator - the initial state is stationary though out-of- 
equilibrium, equilibrium meant to be the two leads at the 
same chemical potential. 

Suddenly the coupling to the bridging region is 
switched on - namely the Hamiltonian changes from 
T + Hint into H = T + Hint + V - and a current starts 
to flow. If 



U(t) 



-mt 



(6) 



is the time evolution operator with the full interaction, 
the initial density matrix po evolves in time maintaining 
the functional form of a Boltzmann exponential 



P(t) 



/Tr(e 



where 



and 



H(t, $) = H(t) + *y(t) 



(7) 



(8) 



H(t) = U(t)(f + H m t)U{t)\ 

Y{t) = U(t)Y U(tf (9) 

For time t sufficiently large, namely after a transient 
time T, the system reaches a steady state with constant 
current. If we are interested only in steady state prop- 
erties, a good starting point is offered by Hershficld's 
results*^ He showed that the stationary state value of 
certain obscrvablcs coincide with their equilibrium value 
obtained through the effective density matrix 



P 



/Tr(e 



with 



H{<f>) = H + &Y, 



(10) 



(11) 



where Y is the asymptotic time evolution of Yq still 
satisfying^ 



[H,Y] = 0. 



(12) 



III. THE RESONANT-MODEL OUT OF 
EQUILIBRIUM 

Let us consider the simple case of a non-interacting 
single-level quantum dot 



H = T + V , 



(13) 



with 



t= E ^2 e * c lk„ c a ko +E e d4^. 



a= — 1,1 ko 



a= — 1.1 ka v 



(14) 



where c aka creates a conduction electron on the left 
(a = —1) or right (a — 1) lead with quantum number k 
and spin a while d\ creates an electron in the dot with 
spin tr, and £1 is the quantization volume of the system. 
Notice that, quite generally only a single channel of con- 
duction electrons is coupled to the impurity, so that the 
model can always be mapped onto two one-dimensional 
leads hybridized at the contiguous edges with an impu- 
rity. Therefore it is perfectly legitimate to regard the 
quantum number k as one-dimensional momentum and 
S7 as the linear size of the system. 

Let us assume that our system does not have bound 
states. In this case it can be prove n 19 i 20 that the nonequi- 
librium Hamiltonian 

Wo($) = Ho + (15) 
can be expressed as 

a— — 1,1 k<7 



a— — 1,1 ka 



where i$ aka are the fcrmionic creation operators that gen- 
erate the left (a = —1) and right (a = 1) incident scat- 
tering waves 



€ka 10) 



1 



e k -U + iO+ 



V 



4. |0) 



E 



~rrz — : , ,- n + (1 7 ) 



being gd(c) the retarded Green's function of the impurity 
at equilibrium, which, in the infinite bandwidth limit, is 
given by 



Should Y be known, steady state properties could in prin- 
ciple be obtained by any equilibrium technique. 



(18) 
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We underline that Eq. (jTT)) is meaningful only in the 
thermodynamic limit, i.e. when D, — > oo. For a finite 
system the time evolution of an incident state 



ICU = 4jo> 



(19) 



oscillates, namely it doesn't converge to a well defined 
scattering state 



loafer) = V»Lal°> ■ 



(20) 



The scattering states (fTT)) constitute, in the thermody- 
namic limit, a complete basis 

E ^akc ^ a k* = E C l^ Qfca + E d ^ d o , (21) 



provided that there exist no bound states^ Eq. ([2"TT) al- 
lows us to formally expand the c and d operators as fol- 
lows 



c -E-=^-f-+E 

aAicr r a/ccr Z > 



a A; 



n e k -e k - i0+ Vaks 



(22) 



and to calculate the average of any operator using the 
result 

^a'W !*(*)> = f (e k + $|) (23) 

- where is the ground state of Ho(^) and /(e) is 

the Fermi function. The correct value of the average is 
finally obtained taking the limit for f2 — > oo of the result. 
It can be proven that the obtained value is the same that 
one could obtain within the Keldish technique. 

It is very important to underline that the scattering 
operators can formally be defined even in the interacting 
case, although their explicit calculation is not feasible in 
practice. If, for instance, we add to the resonant model 
Hamiltonian (fT3|) a Hubbard repulsion term on the im- 
purity 



(24) 



(Anderson impurity model), the scattering operators 
^aktr are defined as the asymptotic time evolution of the 



°afccr 

nian 



operators generated by the full Anderson Hamilto- 



n = T + V + U, 



(25) 



and still satisfy the completeness relation (|2"Tj) in the ab- 
sence of bound states.— Moreover the interacting Hamil- 
tonian (f2"5"j) can still be expressed in terms of scattering 
states^ 



H = 



E 

a= — 1,1 lea 



y= E Eo^T^ 



a= — 1,1 fecr 



(26) 



We conclude this section calculating the energy £ $ for 
the non-interacting model (|13[) . More precisely, we con- 
sider 

(27) 

where £ $ is the energy of the uncorrelated system T. 
Using Eq. (|2"2"j) it can be proven that 



6£* 



' ie n +a— 



ed - A I ie n +a— \ \ 



de 



6(e) 



f e 



/ e 



where e„ are Matsubara frequencies, 



is the hybridization function, and 

'e + i0 + -e d -A(e + i0 + ) 



<y( e ) 



-Im In 



(28) 



(29) 



(30) 



e + i0+ - e d 

If we assume that the half-bandwidth W is the unit of 
energy, that the density of states is flat 

A(*)=/*iM 
J 7T z - e 

r( c ) = r X[ -i,i](e) 

_ / 1 Ve € [-1,1] 

Vc0[-l,l] ' 

and that r -C W = 1, it can be easily verified that 



X[-i,i](e) = | I 



(31) 



2 d_ (5£*{T) 
dT 



2 /r 

- — e arctan — 

7T \ { 



(32) 



We observe that when $ = the right member of Eq. 
is 2F/7T, so that the solution of Eq. (f3"2"j) is 

<5£ = --riog(i 

(being e the Nepero's number), which derives from an 
hybridization gain 



(33) 



5£hyb = 

and a bath energy cost 



7T VI 



2 / 1 

SSbath = - r log — 
7r \ el 



In the regime 

iu = i>$>r, 

the right member of Eq. (f3"2"|) vanishes, so that 
2 

<5£* = -Tlos 



(34) 

(35) 
(36) 
(37) 
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IV. THE BIAS AS AN "ORDER PARAMETER" 

In this section we propose a point of view of the 
nonequilibrium problem in quantum dots based on ideas 
and definitions very similar to those encountered in the 
general theoretical description of collective phenomena in 
quantum mechanics ^ 

In the standard quantum theory of finite systems there 
is a one-to-one correspondence between the observables 
and the operators in a certain Hilbert space. This corre- 
spondence is unique (Von Neumann 1955), and this en- 
sures that the choice of a specific representation of a finite 
system does not lead to loss of generality. The situation 
is different for infinite systems, because Von Neumann's 
theorem can no longer be applied. The observable of an 
infinite systems generally admit a big variety of incquiv- 
alent representations, corresponding to macroscopically 
different classes of states. 

Let us consider, for instance, the case of the Heisenberg 
model 



H — — J S rl S„+i , J > . 



(38) 



A ground state |^) of the system has all the spins aligned 
parallel in the same direction. The class of all the states 
obtained applying a finite number of local spin transfor- 
mations to |$) does not change its magnetic order pa- 
rameter 



m 



lim 



£ (*[s»hg) 

i=-JV 



2N+1 



(39) 



because the system is infinite. Furthermore the topologi- 
cal closure of the space spanned by these states (that is a 
Hilbert space) is the basis of an irreducible representation 
of the algebra of the Pauli spins {S„}. In this example 
we have a set of different phases corresponding to differ- 
ent (inequivalent) representations of the observables for 
any arbitrary direction m. In other words, a state \^ m ) 
defines a corresponding island of states T(^> m ) that share 
the same order parameter and are the basis of a specific 
representation of the algebra of the observables. 

The general formal definition of the islands is based 
on the fundamental Gelfand-Naimark-Segal (GNS) theo- 
rem^ that tells us that if we have an algebra A gener- 
ated by the local observables of an infinite system - i.e. 
all the operators that belong to any finite region of space 
A, - for each state p there is an operator-representation 
A p of A on a Hilbert space Hp, which is determined (up 
to unitary equivalence) by the conditions that it exists a 



vector \^> p ) in Tip such that 



P(A) 



(%\A\%) VAeA, 



(40) 



and that "H p is generated by applying the elements of A p 
to |* p ). 

The island I(p) is, by definition, the set of all the states 
p' corresponding to all the density matrices p' in Hp. 



Physically, the meaning of the states p' £ I(p) is that 
these are states generated by localized modifications of 
p, but are "macroscopically" equivalent to one other. 

Let us now consider the system represented in Fig. [TJ 
whose dynamics is defined by the interacting Anderson 
impurity model 



aka aka 

u 



J2e d dld a + —(h d -l) 



where 



rid 



£44, 



(41) 



(42) 



is the impurity number operator. 

To the initial state po($, T) defined in Eq. (|2|) will cor- 
respond, through the GNS theorem, an island of states 



I($,T)=X(po(*,T)) , 



(43) 



and the islands obtained from initial states with different 
$ will correspond to inequivalent phases, because of the 
infinite volume of the two leads. In this sense we can say 
that $ plays the same role as the magnetic order param- 
eter in the example of the Heisenberg model considered 
above. 

In this work we assume that the dynamics of the sys- 
tem does not mix vectors belonging to different phases, 
i.e. that if p' belongs to Z($,T) then so does p'(t). This 
corresponds to the physical idea that the system never 
equilibrates because the two leads are infinite and the 
junction between them (the dot) is finite, so that the 
current through the dot can not change the densities of 
the two leads defined by the value of $. 

The stability of the dynamics in I($,T) allows us to 
consider the generator of the time evolution trans- 

formation inI($,T), i.e. the Hamiltonian that describes 
the dynamics of the island. 

The basis of the ideas proposed in this paper is that an 
equilibrium problem corresponds to the study of the 3> = 
phase, while a nonequilibrium problem is equivalent to 
study a finite-^ phase. Once the "H^t operator is defined 
we can apply, in principle, any equilibrium technique to 
study the physics of the corresponding island. 



V. VARIATIONAL APPROACH AT T = 

In this section we will concentrate our attention to a 
general GNS representation of the system at T = 



J($) =2T($,T = 0). 



(44) 



In this case we know that the initial state is represented 
by a pure vector 

|* (*)> G X(*) 
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Po($,T = 0) - |* (*))(*o(*)|, 



(45) 



and that the operator that governs the dynamics of X($) 
is simply given by 



1~L^ — 7^$ ; T— o 



(46) 



where P$ is the projector on X($). 

We define Z(4>) as the product of the Hubert space of 
the dot and the Hilbert space generated by |\l/o($)) and 
all the states 

i*05)= n 4jo>, ( 4? ) 

where S is any subset of 

E = {(a, k,a)\a = ±l,k£ [-it, it], a = ±1/2} (48) 
that differs from the set 5$ of the occupied states of 

l*o($)>^ II c ^l°>- ( 49 ) 

(a,fe,ir)(:S$ 

by an arbitrary, but finite, number of particle-hole trans- 
formations. In other words, any state l^os) can be writ- 
ten in the form 

hM= n n c ik^ k -j^m) , (50) 

(a,/c,tr)G/Cs {a. 1 ,k' (7')€lK' s 

where ICs and IC' S are finite subsets of E. 

We are going to prove that Z(<3>) coincides with Z($). 
To prove that X($) C Z(<I>) we observe that whether a 
state |*) £ J($) belongs to also the state 



I* 



aficr.a' k' cr 



E 

RR'eA 



ikR -ik'R' „i 



< Ra c a , R , G ,m (5i) 



belongs to it for any finite region of space A. But, by def- 
inition, !($) is the topological closure of the space gen- 
erated by the states obtained modifying |\l/o(3>)) locally, 
so that the state 



c L CT <W CT ,|v]/) =lim|^ feCT)a , fev ,>, 



(52) 



where liniAf denotes the limit for the size |A| of A going to 
infinity, belongs to too. The inclusion I($) C 
is then proven by induction. 

To prove that !($) C Z(<I>) we consider again the set E 
defined in Eq. (jig]) of all the possible values of (a, k, a), 
and we associate to each subset e C E the operator 



^ A = E 

i a - k =fx\ n - 



. a (n integer ] 



(53) 



Let us consider a general state \ty) £ If |^) does 

not belong to !($) then, by definition, one can define a 
set e C E of single particle states such that either 



or 



lim--(* o (a>)|iV A |vI/ o ( ( f>))=0 
At |A| 



Hmp|(*o($)|^e A |*o(^)>>0 



(54) 



(55) 



But, if such set exists, it is clear that |^) can't be 
generated by the application of local observables to 
l^o because the contribution of local modifications 
to Eqs. (|54j|55[) vanishes in the limit of |A| — s- oo. The 
inclusion I($) C is then proven. 

We can reformulate the statement X($) = X($) saying 
that the measure d/i^(a,k,a) such that 



/ (Wa, k, a) = lim t-tt^I I*) Vec£ 

Je A t |A| 



(56) 



is an order parameter that identifies the phase Z($); i.e. 
that the states I'J) that belong to are characterized 
by the condition 



d/i*(a, k, a) = d/^ ($)(a, k, a) . 



(57) 



We underline the strong analogy between the measure 
d^(a,k,a) and the magnetic order parameter m, see 
Eq. |39|) . for the Heisenbcrg model. 

As we have anticipated in the previous section, we as- 
sume that X($) is stable respect to the dynamics induced 
by the dot. While such assumption is very reasonable for 
finite time evolutions, it is less trivial that the steady 
state, which is reached only after an infinite time, still 
belongs to I( < f>) - and this is what we need. Although 
we can't prove the stability of the respect to the 

asymptotic dynamic induced by an interacting dot, it is 
encouraging to observe that when U = 

e?M*(*) {ct, k, a) = d/i* ($) (a, k, a) . (58) 
This can be easily verified from the following equation 
|*(*))= J] ^ fe j0), (59) 

(a,k,a)£S<i> 

where the scattering operators ip aka are given by Eq. (|17p 
and 5$ C E is defined bv Ea (|4"9"|). 

As a consequence of the stability of the (asymptotic) 
dynamics induced by the dot we can characterize I( < I ) ), 
in the absence of bound states (when Eq. (T2T|) is satis- 
fied)^ even as the space generated by | V I'( < I > )) and all 
the asymptotic time evolutions of the eigenstates of T-Lq 
defined in Eq. (|4"T)) . i.e. the states 



i*s>= n ^alO), 

(a,fe,cr)eS 



(60) 



where S differs from the set iS$ of the occupied states of 
|^($)) by an arbitrary, but finite, number of differences, 
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and il^ ak(J are the scattering operators of Eq. ([25)1 (that 
are interacting in general). 

Starting from the above characterization of we 
can understand that all the possible eigenstates Y&s) oiii 
defined in - i.e. the eigenstates of - correspond 
to the same current 



Is = -i 



ka 



^<*($)l4 Cfca> _ 1 1*(*)>-< 



(61) 



and can be considered "equivalent" in this sense. Any of 
the eigenstates \H> S ) of H in T($) defined in Eq. (JHOJ)- is, 
in fact, the asymptotic time evolution of the correspond- 
ing "initial" state 



I* 



OS) 



n 



(a,k,a)eS 

The existence of a particular set S such that 



(62) 



(63) 



would imply that the state \^qs) ~ obtained, by defi- 
nition, modifying the initial state |\&o($)) applying to 
it only a finite number of particle-hole transformation - 
leads, after an infinite transient time, to a different cur- 
rent respect to the one of the steady state |^($)), and 
this is clearly unphysical (notice that the same argument 
can be applied to any local observable, and not only to 
the current operator). 

It is interesting to check directly the validity of the 
above statement for the simple non-interacting case U = 
0. The scattering operators are, in this case, given by 
Eq. (jTT)) . and the average of the current operator can be 
calculated with Eq. (|2"2"j) and the identity 



where 



|*s)= <W<W<W n^ s (a, k, a) ; (64) 



n$ s (a, k, a) = 
The result is 



1 if (a, k, a) € S[#] 
if (a, k, a) <£ S[*j 



(65) 



(66) 



(67) 



(68) 



i.e. the difference between the two functions does not 
contribute to the integral that defines the value of 
the current /, because, by definition, 



1=1 der(e)p d (e)^2an^ s (a,e,a) , 

^ eta 

where pd(e) is the spectral function of the dot 

Pd(e) = --lm(g d {e + i0 + )) . 
If |*) el(l) = !($) then 



n* s (a,e,cr) ~ / + $ 



E 

Q>(7 ' 



de 



nn, s (a,e,a) - / ( e + $ 



= 0. 



(69) 



In order to identify the eigenstates of ji® a possibility 
is to minimize the variance 



4 [n?] 



(<f (<s>)\n^Hlnf\* (<t>)) 
(HF {$)\n<?si°\v ($)) 

( <^ (^)| ng^ng |^o(#)> ' 
1 (* (*)|nfn?|*o(*)) 



(70) 



respect to the most general operator f2? that is gener- 
ated by the algebra of the local observables. Notice that 
what we have defined is a variational principle for the 
Hcrshfield steady state at zero temperature! 



A. An energy-based approach 

One may be tempted to use the energy instead of the 
variance, i.e. to claim that the ground state can be cal- 
culated variationally even minimizing the energy 



(* (*)|n+ t n?l*o(*)> 



(71) 



respect to the fl ^-operators defined above. 

Unfortunately such energy minimum does not exist, 
because a $ > phase contains states obtained from the 
Hershfield state | V I , ( < I ) )) moving an arbitrary (although 
finite) number of electrons from one of the leads to the 
other. A solution to this problem would be to consider 
only the limited subset of f2+ such that fi+|\I>o($)) satisfy 
the equation 



^ng|* (*)> = flV ($)fi?|¥ (*)> 



where 



Y 



E 

a— — 1,1 ka 



Ef ^ 



aha 



and SNq($) is defined by 



f |* ($)) = <Wo(*) !*„(*)) 



where 



ka ' 



(72) 
(73) 

(74) 
(75) 



a= — 1,1 ka 



In this subspace it is clear, if we think in terms of in- 
teracting scattering operators, that the only eigenstate 
of "H$ is the state that minimizes the energy. Unfortu- 
nately the operator Y is not known, so that the energy- 
based procedure defined here can't be applied rigorously 
in practice. 

Notice that if we formally apply the exact (unknown) 
asymptotic time evolution operator f2+ to the two mem- 
bers of Eq. ([71)1 



fi + y |*o($)) = w ($)o + |$o(f)) 



(76) 
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we obtain exactly the condition defined in Eq. (1721) . 
namely that 

ra+ 1 *„(*)) = (n+ nt) n+|*o(*)) 

= 5JV (*)n+|*o(*)>, (77) 
where we have used that 

Y = 0+ Y a nl . (78) 

If we apply the approximated trial f2+ to Eq. ([71)1 we 
obtain, instead, that 

Slf f |* ($)) = <WV ($) fi£ |* ($)) , (79) 

which is equivalent to 

y g |* (*)} = <Wo($) |* (*)} , (so) 

where 

Y a = n%Y a°- 1 . (8i) 

For this reason, whether we believe, for some physical 
reason, that the proposed fi9 is sufficiently good to guar- 
antee (approximately) the equivalence of Eq. (|80[) and 
Eq. (|72p the energy minimization procedure is still mean- 
ingful, although not purely variational. 

In the following sections we will propose an example 
of such energy-based approach starting from a particu- 
lar variational space, that is expected to describe suf- 
ficiently well the qualitative behaviour, in the Fermi- 
Liquid regime, of the interacting single-orbital Anderson 
impurity model. 

We conclude this section observing that, at least for- 
mally, the GNS theorem mentioned above allows us to 
define a variational principle for the Hershfield steady 
state even at finite temperature. In fact, to a thermal 
state po(&,T) corresponds, in its GNS representation, a 
vector \tyo($>,T)), and the dynamics of the correspond- 
ing island X($, T) is (presumably) still generated by some 
operator H<s>,t, although it is not simply a "block" of H 
(Eq. ((46])). An eigenstate of H<s> y T is a stationary state, 
i.e. the representation of the Hershfield state in X($, T), 
and can be formally identified by the minimum- variance 
condition as before. 



VI. THE CONCEPT OF QUASI-PARTICLES 
OUT OF EQUILIBRIUM 

Let us consider again the general interacting system 
described by the Hamiltonian 

H = Ho + V + n in t, (82) 

We know that if we prepare the two leads at a different 
chemical potentials and let it evolve within the interact- 
ing Hamiltonian 

U{t) = e-***, (83) 



for times t longer than some transient time T the fi- 
nal noncquilibrium state is described by the Hershfield 
Hamiltonian 

H($)=H + $Y, (84) 

formally defined in Eq. 

In general Y is a complicated many body operator that 
must satisfy Eq. (|12j) and in addition share the same sym- 
metry properties as Yq, i.e. a spin-singlet operator odd 
under interchanging the two leads. Therefore, generally 
the steady-state Hamiltonian %($) is an interacting one, 
the interaction 

5H($) = Hint + * (y - %) (85) 

presumably remaining "local" (in the sense defined in 
section |V} as it was originally. Furthermore, since the 
nanocontact can not change the bulk properties of the 
leads, e.g. inducing a spontaneous symmetry break- 
ing, H(&) should still describe a metal. It is there- 
fore tempting to assume that, if in the absence of ex- 
ternal bias the system, leads plus nanocontact, is de- 
scribed by a local Fermi liquid theory in the Nozicrcs 
sensed which is generally the case, the same should 
hold even in the steady state after the bias is applied^ 
It then follows that it should be possible to represent 
the low energy/temperature/bias properties in terms of 
weakly interacting quasi-particles which, by continuity 
with the non-interacting case, should be better regarded 
as renormalized scattering states with an Hamiltonian 
of the same form as (|I6[) with renormalized (bias de- 
pendent) energies plus additional weak local-interaction 
terms' This local Fermi-liquid assumption seems to us 
quite plausible. However, since the bias is coupled to 
a non-conserved quantity, the charge difference between 
the leads, the effective bias felt by the quasi-particles 
will generally differ from the applied one and the quasi- 
particle current does not correspond to the real one. This 
implies that the current can not be expressed simply in 
terms of Landau parameters and an explicit calculation 
is required. 

VII. THE GUTZWILLER APPROXIMATION AT 
EQUILIBRIUM 

Let us consider the 4> = phase, namely the equilib- 
rium problem. Although the method we shall present is 
quite general, for sake of simplicity we shall show how it 
works in the simple case of a bridging region described by 
a single-orbital Anderson impurity model at half-filling 




+ -{n d -lf = (f + y) +11 = Tio + U (86) 
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FIG. 1: (Color online) The single dot system. 



The physical properties of the above Anderson impu- 
rity model are very well known^ For large U the model 
effectively maps into a Kondo model, the impurity elec- 
tron behaving as a local moment Kondo screened by the 
conduction electrons. A simple way to describe qualita- 
tively and to some extent also quantitatively the Kondo 
screening is by a Gutzwiller-type of variational wavefunc- 
tion^^ 

|*>=P d |*o) (87) 

where Vd is an operator that modifies the relative weights 
of the impurity electronic configurations with respect to 
the uncorrelated wavefunction \^o) 

Vd = A 0) 2 (|0)(0| + | U)(U I) 



+ Ai(|tXt I + UX^|) , (88) 

and l^o) is the ground state of a non-interacting varia- 
tional resonant level Hamiltonian. 



A. Energy optimization 

The variational procedure amounts to optimize both 
the local projector Vd as well as the non-interacting wave- 
function l^o) by minimizing the expectation value of the 
Hamiltonian 

We assume that Vd is subject to the following two con- 
ditions 

(*o|P]P d |*o) = I. (89) 
(^o\VlV d h d a\^o) = (*o|*W|* >, (90) 

where 

hda = d] J d a . (91) 

Condition (|89|) is the normalization requirement of the 
variational wavefunction, that corresponds, in terms of 
A-paramcters, to the condition 

^o,2 + A? = 2 . (92) 
Condition (|9"U)) - that ensures that all the Wick contrac- 
tions between the conduction electron operators and the 
impurity operators are zero - allows to evaluate expecta- 
tion values straightforwardly. 

In particular, the expectation value of the Hamiltonian 
(|8"6")l . that has to be minimized, is 



(V\H\V) 



cxka aka 



l*o> + ^(*o|^(^-l) 2 ^|*o> 



(93) 



where the hopping rcnormalization coefficient R is ob- 
tained through the following equation: 

(*o|^4^^l*o> =i?(* |4d ff |*o), (94) 
whose solution is 

R = A , 2 Ai . (95) 

The calculation of the first term in Eq. ((93)) reduces, pro- 
vided eqs. (|89[) and (|90l) arc satisfied, to calculate the 
energy gain of 7-^ due to the renormalized tunnelling 



term 

V* = E^I^W+^c (96) 

The variational Hamiltonian whose ground state is the 
uncorrelated wavefunction |Wo) has rigorously no phys- 
ical meaning but for the ground state properties. How- 
ever, it is common — to interpret it as the Hamiltonian 
of the quasi-particles and 

R 2 = z (97) 
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FIG. 2: (Color online) z, as a function of U/Y. 



as the quasi-particle weight of a single-particle excitation. 
Within such an assumption, the Gutzwiller approxima- 
tion technique can be regarded as a tool to extract quasi- 
particle properties. 

From now on the unit of energy is given by the con- 
duction electron half-bandwidth W . The explicit value 
of the variational energy (referred to the ground state 
energy of the unperturbed system), is then given by 

6E[V] = ^zT log (^j + j (1 - VT^) ■ (98) 

In Fig. [2] we show the value of the optimal z, as a 
function of U. At U — wc find that z — 1 (as expected) 
and has a finite curvature. When U — > oo we find that 

z (c/)~Iex P (-^r9 • (99) 

Notice that at large U the value of z vanishes expo- 
nentially but remains finite because 

(*o|#§il*o> = --*riog(4) , (ioo) 

7T V Zl / 

which vanishes at z = with an infinite derivative due 
to the presence of z in the logarithm. 

We conclude this section by underlining a limit of the 
Gutzwiller method from a quantitative point of view. Us- 
ing Eq. we can define the "Gutzwiller approxima- 
tion" for the Kondo temperature as 

^~^)^exp(-^). (101) 

Notice that Tjf differs with respect to the correct value 
of the Kondo temperature 

TK-expt—tpj (102) 

because: 



• the universal prcfactor in the exponent should be 
7r/8 and not 7r/16, 

• the factor W/T (which is equal to 1 /T in our units) 
in Eq. t|101[) diverges in the infinite bandwidth 
limit. 

The divergence of the right member of Eq. (|101[) for 
W/T — > oo reflects the unreliability of the method in 
this limit. In fact, when W/T — ¥ oo, the Gutzwiller ap- 
proximation predicts that z — ¥ 1 even if U > F, as can 
be verified directly from Eq. (|9"5j) . 

B. Variance optimization 

For the equilibrium ground state of a system it is pos- 
sible to define an infinite number of functionals with the 
same minimum. For example, the functionals 

F^]=E[^]+ f ia H m, (103) 
being a-u the variance 

**M = [Mw 3 M-M«h&> 2 ]' , (104) 

are equivalent V/x > 0. Nevertheless the result obtained 
when the functional is minimized on a particular vari- 
ational space depends on the factor /z, and the choice has 
to be motivated on the basis of the specific problem con- 
sidered. 

In particular, if one think that his variational function 
is a good approximation of the ground state (and not of 
the excited states) of the system, instead of minimizing 
the energy E, it is sometime convenient to minimize the 
variance cr-^. From the equilibrium variational principle, 
the smaller the energy is the better the variational state 
will be, but, without an exact solution, it is hard to judge 
how accurate the variational approximation is. On the 
contrary, the variance is very useful, because the small- 
est possible variance, equal to zero, is known a priori, 
and in this case the variational state represent an exact 
eigenstate of the Hamiltonian. 

We observe that in the Gutzwiller energy-minimization 
procedure discussed in the previous section the optimiza- 
tion scheme leads to a correspondence between the A- 
parameters (defining Vd through Eq. ([55])) and the uncor- 
rected wavefunction l^o): once the A-parameters were 
defined the corresponding |^> ) was the ground state of 
a renormalizcd "variational Hamiltonian" "H^j, with R 
given by Eq. ([55)1 . In the appendix we show the calcula- 
tion of the variance assuming such correspondence. 

If, for simplicity, we assume that 

V k = V k > Vfc, k' 

r < w = i , (105) 

the variance is given by. 

= ( i-,)g + ^r 2 ziog 2 (z 2 r 2 )) 
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u 



VT^lr log (z 2 r 2N > 



zU 2 



(106) 



Unfortunately, the minimization of the variance func- 
tional (|106[) respect to the allowed values of z does not 
lead to a physically reasonable result (not shown). Our 
conclusion is that the minimization of the variance re- 
quires, to be effective, a more realistic trial state respect 
to the simple Gutzwiller-type wavefunction (|87[) consid- 
ered here. 

Following the derivation of Eq. (|106|) in the appendix it 
is clear that only for Gutzwiller's wavefunction such that 
the correspondence ([§5]) between |\I>o) and Vd is verified 
the variance is finite. We underline that at half-filling 
such correspondence eliminates the diverging terms (for 
our system) even out of equilibrium, i.e. when |$o) is 
the ground state of 

aka 
aka 



instead of the ground state of 



Hi = Y efc c f , c , + V ^ 4c , + H.c. 

tt / j k aka aka / j akcr 
akcr akcr 

= E £ *^( fl )^(fl); ( 108 ) 



akcr 



being 



RV, 



k R 

9d 



(e*)4 



ot'k'o 



^ R?V k V k , g«{e k ) t 



the appropriate rcnormalized scattering operators iden- 
tified by R. 

This observation suggests that at half- filling the simple 
form of the trial function (|57| is a reasonable variational 
representation not only for the ground state of our sys- 
tem, but also for its noncquilibrium Hcrshfield steady 
states. This observation relates with the Fermi-liquid as- 
sumption formulated in Sec lVIl 



VIII. THE GUTZWILLER APPROXIMATION 
OUT OF EQUILIBRIUM 

We study now the half-filled Anderson model 



H = E £fc cl ka c aka + E -7= d U ak a + H 



aka 



akcr 



(110) 



when it is driven out of equilibrium preparing the leads 
at two different chemical potentials (Fig. [1| . Turning on 



the tunnelling interaction we know that a current starts 
to flow and the system, after a transient time, reaches 
the steady state formally defined by Eq. (jTOJ) . At zero 
temperature the steady state is therefore the ground state 
of 



U(<f>) = H + $Y. 



(Ill) 



We want to approximate the Hcrshfield steady state 
with the usual equilibrium Gutzwiller variational wave- 
function 



I*) =^1*0) 



(112) 



which satisfies conditions and (|i?0"|) . 

The average on |\I>) of the non-interacting part of H 

Wo = E e * c lk^ aka + E "If 4c ato + H.c. (113) 

aka aka 

is equal to the average on | ^0) of the rcnormalized non- 
interacting Hamiltonian 



o}0 



aka 



E^4w+^.c(ii4) 



The difference between the equilibrium and the noncqui- 
librium case is that in the presence of a bias we can con- 
sider only Slater determinants belonging to In 
order to guarantee the (approximated) equivalence of 
Eq. (|5U|) and Eq. (fT2"j) the only reasonable state is the 
Hershfield steady state of the renormalized uncorrelated 
system (|114[) . namely the ground state of the Hamilto- 



(115) 



defined in Eq. (fTUTJ). In fact, just like at equilibrium, 
we expect that the Hamiltonian (|107p should describe 
weakly interacting quasiparticles in the presence of a bias. 
Our particular choice for the variational Slater determi- 
nant l^o) hi Eq- (|H2|) simply means that the number of 
left and right quasiparticle is equal to the number of left 
and right particles in the unperturbed system. 

Summarizing, our variational choice corresponds to ap- 
proximate the asymptotic time evolution operator f2+ 
with 



^+ — Vd Q+r , 



(116) 



being OVr the unitary operator (to be determined vari- 
ationally) such that 



(117) 



The corresponding procedure amounts to minimize the 
following energy functional 

= (* (*)|#hI*o($)} 
+ % (*o(*)| Vl(h d - l) 2 V d |* (*)> , (118) 
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FIG. 3: (Color online) z as a function of U/T for V 
and three different values of the bias V/T. 
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5x10" 



where |* ( < i ) )) is the ground state of Ti R (^) (that sat- 
isfy the conditions 



(|89H90p ). In other words, the only 
difference respect to the equilibrium calculation is that 
we substitute the equilibrium energy gain due to the tun- 
nelling term (|96|) for the energy gain due to the tunnelling 
term in the non equilibrium quasi-particle Hamiltonian 
(|107p . It can be easily proven that the value of the vari- 
ational energy (referred to the noncquilibrium energy of 
the unperturbed system), is given by 



U 



6E*[V] = + j (1 - VI 



(119) 



where 



2 . ^ N $ 

■— arctan(zr) H arctan . 

7T 7T V^/ 2 



(z^ 



— log 

7T 



z 2 r 2 



($/2) : 



i + z 2 r 2 



(120) 



We stress that our functional, and then the value of R 
after the optimization, depends on the bias 4>. This is 
crucial in order to properly take into account the strong 
correlation effects induced by the Hubbard repulsion, i.e. 
to obtain the expected destruction of the Kondo reso- 
nance at finite bias (see Fig. [3]). 

The expression for the average of the current after the 
optimization is 

1 = -^E"^ 4c fc(rt _il*o(*)>-c.c. 



fe<7 



der R (e)p r d *(e) 



(121) 



where p d R (e) is the spectral function of the dot, that is 
1 T R 



r 



X[-i,i]( e ) 



(122) 



with 



T R (e) = R 2 T(e) 



(123) 



having assumed that the density of states is flat and that 
r < W = 1 



J tt z-e 

r(e) = r X[ _ lil] ( e ) 



(124) 



We notice that Eq. (|121[) fails to describe the sys- 
tem accurately when $ ~ U, because it doesn't take 
into account the spectral contribution of the Hubbard 
bands. However, for the simple single-band Anderson 
model we can reproduce artificially the correct qualita- 
tive behaviour of the current in this regime by substitut- 



ing R 2 p T d R (t) with 



a=-l,l 



e-U- 



(125) 



in Eq. ([121)1 ■ 

In Fig. [D we show the results for the conductance G 
of the Anderson model. The obtained value of the con- 
ductance at zero bias is universal as expected, and the 
curvature is given by 



d 2 G 
d<P 2 



1 



1 



27r(i? 2 F) 2 (T§ ) 2 



(126) 



- Tj£ being the Kondo temperature with the incorrect 
prefactor predicted by the Gutzwillcr method 



T, 



K 



(127) 



Nevertheless for large enough value of U we found (not 
shown) that the conductance may become negative, 
which is unrealistic. In order to establish the regime of 
validity of our method, we note that the Fermi-liquid de- 
scription that we assume is applicable only for values of 
the bias much lower then the Kondo temperature Tk- 
For the single-orbital Anderson impurity model we can 
calculate analytically the minimum value of the energy 
functional (|118[) when 



namely when Eq. Q37[l can be applied, so that 



(128) 



(*o(3>)|^ + ^|*o($)> = ^i? 2 riog^|^ . (129) 

In particular, it can be easily proven that the value of z 
vanishes at 



~2 



e 16 r 



Ik, 



(130) 



that is out of the expected regime of validity of the cal- 
culation. 



12 





u/r=2 




u/r=6 




u/r=io 



2 4 6 8 10 12 14 16 18 20 

v/r 



FIG. 4: (Color online) Conductance as a function of the bias 
V/r for T = 5x 10" 3 and three different values of U/F. 



A. Why is half-filling special? 



We observe that, if /i is taken to be the value that satis- 
fies the constraint (|^TJ|) at equilibrium, the ground state 
(^)) °f the non-equilibrium Hamiltonian 

aha 

+ *Ef^( fl )^w ( 137 ) 

afccr 

is not such that 

(tt (*)l44rl*o(*)>=», (138) 

namely it doesn't satisfy anymore (1901) . The procedure 
described for e d = should then be modified without 
fulfilling this condition. This forces us to renounce to 
the mechanism that eliminates the diverging terms of the 
variance mentioned in Scc. lVII TH and that we interpreted 
as measure of merit of the variational ansatz defined in 
Eq. (fTTol) . 



Let us consider the Anderson model away from 
particle-hole symmetry 



akcr a.ka 



e d ^2dld a + — {n d - 1) 



(131) 



by requiring e d ^ 0. The state |^o) which minimize the 
energy n° R and satisfies Eq. (|90j) 

(* |4dJ*o>=n (132) 

can be calculated within the Lagrange multipliers 
method, namely is the ground state of the Hamilto- 
nian 



(133) 



with a proper chemical potential \i. 

In particular, when e d = the ground state of TL R 
satisfies the constraint (|9"0"|) automatically, namely H R = 
"H R , and the corresponding non-equilibrium Hamiltonian 
H R ($) automatically satisfies the constraint (|90| too, 

(tt„| d\d a |tt ) = (*o(*)| 4d„ |*o(*)) = \ (134) 

Let us now consider the general case e d ^ 0. In this 
case 



n R = ^2e k ^ aka (R)^ aka (R), 



(135) 



aha 



where ipl,u a (R) where V'lfe cr (-R) are the scattering waves 
constructed with renormalized hybridization R Vk , which 
depend on the retarded impurity Green's function 



9dR(e) 



e - /it + iR 2 T 



(136) 



IX. CONCLUSIONS 

We have defined a variational principle based on 
the minimization of a bias-dependent functional of the 
variance for studying the steady-state zero-temperature 
properties of a general quantum-dot driven out of equi- 
librium through the application of a bias. We have pro- 
posed a similar (although approximated) method based 
on a "constrained" minimization of the energy. 

The ideas proposed in this paper are mainly inspired 
by the Hcrshfield's point of viewi^ that the out-of- 
cquilibrium steady state can be regarded as the equi- 
librium one with an Hamiltonian %($) that includes an 
effective nonequilibrium term proportional to the bias 
$y. Our main result is that, equivalently, the steady 
state can be identified by the following conditions: 

• the initial state po(&) (identified by $lo) defines a 
phase !(<!>) of "macroscopically equivalent" states 
that contains also the steady state p( ( f ) ) 

• jo($) is the "only" stationary state of the correlated 
Hamiltonian H in !($)• 

Such characterization of the nonequilibrium steady state 
does not require the explicit knowledge of the Hersh- 
field's operator Y, and, for this reason, we believe that it 
constitutes an interesting and useful formulation of the 
problem. It opens, in fact, a new possibility to treat 
nonequilibrium using equilibrium methods. 

In order to test our methods, we have considered the 
simple single orbital Anderson impurity model at half- 
filling, finding a good qualitative accord with the ob- 
served behavior in quantum dots for the expected regime 
of validity. The choice of the variational space was, in 
fact, based on the the assumption that the effective Her- 
shficld Hamiltonian H($) describes a local Fermi liquid 
theory^ in the Nozieres sensed 
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The ideas that we have proposed have the big advan- 
tage of being very simple, and we believe that further 
developments will enable us to deal with more compli- 
cated situations and variational spaces. 
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J 



The variance of the Anderson impurity model (|86[) re- 
spect to our Gutzwillcr variational function (|87[) is given 
by 



(*o\PlT 2 V d \y )-(* \vlTV d \y } 2 
\*o\Vl TVV d |* ) + (*o| V\ VTV d |*„) - 2 (* | V\ TV d |tf )(*o| V\ VV d |tt >) 
'(*o I V\ V 2 V d | *„) - (* | V\ V V d |* ) 2 ) 

"(*o | V\ fti V d \* ) + (* 1 V\ tit V d \* ) - 2 (* 1 V\ t V J* > <*o I V\ UV d \ * )) 
(*o| V\ U 2 V d |* ) - (*o| V\ UV d |* ) 2 ) 

I 



(A.l) 



Our calculation can be considerably simplified by the 
following considerations. 



• The variance of the renormalized Hamiltonian 

K = T + V R (A.2) 
on the function |Wo) 

[*o] = ((* o| T 2 1* ) - (* | f |* > 2 ) 
+ R^o\fV\^o) + (^o\Vf\^ Q ) 

- 2(* |T |*o)(*o|^|*o>) 
+ R 2 { (* | V 2 |* > - (*o| V |*o) 2 ) (A.3) 



is zero, because |\?o) is, by definition, the ground 
state of Ho- 
rn Our variational function satisfies the Gutzwillcr 
constraint defined in Eq. (|90p , so that the following 
equations holds: 



(*o|^yPj*o) = R{9o\V\9 ) 
• A direct calculation shows that 

(^o\VlfVV d \^ ) =R(V \fV\* ) 



(A.4) 



(^o\V d VfV d \^ ) =R(^o\Vf\^ ) (A.5) 

{* \Vl TUV d |* ) + (*o| ^ UTV d |* ) 

-2 (* |P]TP d |* )(*o|P] f>^ d l^o) =0 (A.6) 

(*o| t/*7P d |* ) + (*o| Pj UVV d |* > 



-2(*o|P^P d |*o)(fo|P]FP d |* ) 
AiAo, 2 (1 - Ag i2 ) (* | |*o> 



U 



^ \V^ d U 2 V d \^o) - (^o\VlUV d \^ ) 



u\ 2 a 2 



Aq.2 



(A.7) 



(A.8) 



Taking the difference between Eq. (jA.lj) and Eq. (|A.3j) 
and using the above equations we obtain that 



4 [*] 



(*o|Pl^r 2 |* )-(* |r 2 |* ; 



^ AiAo^ (1 - Ag 2 ) (* |1> |* ) 



c/\ 2 AL 



1 - 



x 2 

A 0,2 



(A.9) 



14 



Let us consider now the first term in Eq. (|A.9[) . which 
is equal to the sum of all the Wick contractions in which 
the operators V^V and T 2 are connected by two or four 
"legs". 

A direct calculation shows that 

V\V d = Ag )2 -(Ag )2 -A;)(4d t + 4d ; ) 

+ 2 (\ 2 02 - A 2 ) d\d^ 4<*4. . (A.10) 



Using Eq. (|A.10[) it can be easily verified that the sum 
of all the two-legs contraction between V][Pd and T 2 is 
zero. The four legs contribution can be calculated using 
Wick's theorem. The result is 

(*o|P]P d T 2 |*o>-<*o|T 2 |tf ) 

= (*h-^)A 2 (V d ), (A.ll) 



1-R Z 



U 



. ka 



R 2 



u\ 2 a§ i2 



1 



Aq,2 



2 J 2 \ 2 



(A.15) 



Notice that the functional relation ([55)1 between |^o) 
and has been responsible of the cancellation of the 
TV terms in Eq. (|A. 1|) . that are extensive quantities. 

It can be easily proven that the minimum condition of 
a^ Q in the not projected state implies that 



where 



A(V d ) = J2^\{^ \4c ak( r\^o)\ ■ (A.12) 

kaa 



Let us consider the second term in Eq. (|A.9j) . It can 
be verified that 



A 2 (V d ) = 0. 



(A.16) 



The explicit value of V(Vd) can be simply obtained re- 
placing T with zT (see Eq. ([57])) in Eq. ([53]). 
If. for simplicity, we assume that 



Wk V V d |* ) - R z <*o| V\ V 2 V d \* 



where 



(A.13) 



(A.14) 



The above calculations lead to the following expression 
for the variance: 



V k =V k , Mk,k' 

r< w = l, 



(A.17) 



it can be easily verified that the variance is given by 



= ( i-,)g + ^ r ^iog 2 (z 2 r 2 )) 



f7 



- z^/^^ r log (z 2 r 2 ) + — , (A.18) 

7T 16 



= (AL-A 2 )^ 2 (P d ) 



that coincides with Eq. (| 10611 . 
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